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The techniques of microwave quantum optics are applied to collective spin excitations in a macro- 
scopic sphere of ferromagnetic insulator. We demonstrate, in the single-magnon limit, strong cou- 
pling between a magnetostatic mode in the sphere and a microwave cavity mode. Moreover, we 
introduce a superconducting qubit in the cavity and couple the qubit with the magnon excitation 
via the virtual photon excitation. We observe the magnon-vacuum-induced Rabi splitting. The 
hybrid quantum System enables generation and characterization of non-classical quantum States of 
magnons. 
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I. INTRODUCTION 


II. MAGNONS IN FERROMAGNET 


The successful development of superconducting qubits 
and related circuits has brought wide opportunities in 
quantum control and measurement in the microwave do- 
main m- In Circuit quantum electrodynamics and mi¬ 
crowave quantum optics, bosonic excitations of the elec- 
tromagnetic modes, i.e., “photons” are handled with high 
accuracy [3{Tn][n]- Therefore, it is natural to extend 
the targets to other quantum mechanical degrees of free- 
dom. The examples are found in recent reports on hybrid 
quantum systems based on superconducting circuits: For 
example, paramagnetic spin ensembles nans] , nanome- 
chanical oscillators [HHE], and surface acoustic waves in 
a piezoelectric substrate m, have been coherently con- 
trolled via a coupling with a superconducting qubit. 

Our goal here is to apply the techniques of microwave 
quantum optics to collective spin excitations in ferro- 
magnet. Similar to superconductivity, ferromagnetism 
has a rigidity in its order parameter. The lowest en- 
ergy excitations are long-wavelength collective spin pre- 
cessions. We couple the quantum of the collective mode, 
a magnon, to a microwave cavity as well as a super¬ 
conducting qubit to reveal its coherent properties in the 
quantum limit [nnn]. 

This paper is struetured as follows: Section 2 reviews 
the basies of magnons in ferromagnet. In Sec. 3, hy- 
bridization of a magnon and a photon in a microwave 
cavity is demonstrated. Finally, in Sec. 4, we demon¬ 
strate strong coupling between a superconducting qubit 
and a magnetostatic mode in a ferromagnetic crystal. 
The magnon vacuum induces Rabi splitting in the qubit 
excitation. Summary and outlook are presented in Sec. 5. 
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A. Spin waves 

In order to describe spin waves, or collective excita¬ 
tions in ferromagnetic materials, we begin with a simple 
Hamiltonian: 

= ( 1 ) 

* ihj) 


where the first term represents the Zeeman energy and 
the second one is the nearest-neighbor exchange interac- 
tion. The sum in the second term is taken over the pairs 
of the neighboring spins. is the Heisenberg spin op¬ 
erator for the i-th site, g is the ^-factor, /iB is the Bohr 
magneton, Bz is the static magnetic field along the 2 : axis, 
and J is the exchange integral. J takes positive values 
for ferromagnetic materials, leading to the ferromagnetic 
ground state, where all the spins are aligned along the ^ 
axis. 

We can express the Heisenberg operators in terms 
of the bosonic operators q, c \ by using the Holstein- 
Primakoff transformation [20] : 
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where s is the total spin on each site. The meaning of 
this transformation is illustrated in Fig. We find from 
Eq. @ that the number of the bosons corresponds to the 
reduction of the z-component of the total spin. 

The bosonic operators defined on each lattice point are 
related to the spin-wave operators by the Fourier trans- 
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which gives the quadratic dispersion relation in the long 
wavelength limit: 

ficjk = 2sJal\k\^ + (10) 

As indicated in the first term, The rigidity of the ordered 
spin System lifts the degener acy of the spin excit at ions, 
which is in stark contrast with the case in paramagnetic 
spin ensembles. 


B. Magnetostatic modes 


FIG. 1. Relation between magnon number and z-component 
of total spin. The number of the bosons corresponds to the 
reduction of the total spin. 
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where N is the number of the atoms with spin 5, and 
Ck and cj^ correspond to annihilation and creation of a 
magnon in the plain-wave mode, respectively. Substitut- 
ing these operators into the Hamiltonian [Eq. Q] and 
truncating it to the second order, we obtain the spin- 
wave Hamiltonian: 


n = yfkJkclck, ( 7 ) 

k 

with the dispersion relation: 

Swk = 2 sZJ (1 - 7k) + QUbBz, (8) 


where Z is the coordination number of each site. In the 
case of a simple cubic lattice {Z = 6) with the lattice 
constant uq, 7k becomes 


7k = 7 (cos kxa^ + cos kya^ + cos kzao), 
o 


While we assumed infinite lattices in Sec. ra we ac- 
tually need to consider finite samples coupling to the 




FIG. 2. Geometries of the transverse magnetization m in a) 
the Kittel [(1,1,0)] mode and b) the (4,3,0) mode. 


surrounding electromagnetic field. Especially when the 
wavelength of magnons is comparable to the sample size, 
the effect of the dipolar field generated by the spins be¬ 
comes dominant, and thus the boundary conditions at 
the sample surface are of great importance. Relatively, 
the contribution of the exchange interactions (oc |kp) be¬ 
comes negligible. We then determine the magnetization 
oscillation modes from classical electrodynamics. 

Suppose the magnetization in the sample is forced 
to oscillate by a time-dependent magnetic field 
perpendicular to the static magnetic field, where 
h = (hx^hy^O). Then the oscillating part = 

(m^,, m^, 0 )e*^^ of the magnetization obeys the Landau- 
Lifshitz equation: 
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where Mg is the saturation magnetization, /io is the vac¬ 
uum permeability, and 7 = is the electron gy- 

romagnetic ratio. Note that we have linearized these 
equations assuming that the amplitudes of m and h are 
small compared to Mg and Bz/ respectively. Solving 
Eqs. ( 11 ) and (12) for rrix and niy yields the following 
results: 
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where is the scalar potential of h = /ioh). Sub- 
stituting Eqs. (13), (14) into the Maxwell equations: 


= /io div h = —div m, 


(15) 


we obtain the differential equation for ip inside the sam- 
ple: 


(1 + k) 
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Outside the sample, = 0. The boundary conditions 
are the continuities of ip and the normal component of 
h + m//io. Walker solved these equations for spheroidal 
samples and found for each eigenmode the mode shape 
and the 5;2-dependence of the eigenfrequency EiHiai. 
These modes, bunched in a frequency range and charac- 
terized by three integer indices, are called magnetostatic 
or the Walker modes. The simplest is the (1,1,0) mode, 
or the so-called Kittel mode, where all the spins in the 
sample precess in phase and with the same amplitude. 
The geometries of the transverse magnetization m in the 
Kittel mode and one of its degenerate modes [(4, 3, 0)] 
are shown in Fig. 


is responsible for the broadening m- The effect of such 
mechanism also diminishes at lower temperat ures. 

A few relaxation mechanisms remain even at the zero- 
temperature limit. For instance, van Vleck’s theory [28] 
assumes an interaction of magnons with an ensemble of 
two-level Systems (TLSs) which have the same excitation 
frequency as magnons. The relaxation rate is predicted 
to be proportional to tanh(fkJni/2/cBT), where /cb is the 
Boltzmann constant. The characteristic temperature de- 
pendence is derived from the fact that the TLSs are sat- 
urated as the temperature increases. To the best of our 
knowledge, however, there has not been any observat ion 
of such temperature dependence in ferromagnetic reso- 
nance linewidth. 

Another relaxation mechanism independent of the 
temperature is the elastic scattering of the Kittel-mode 
magnons due to the surface roughness of the samples [29] . 
The surface roughness causes intermode coupling be- 
tween the Kittel mode and other magnetostatic and spin- 
wave modes. 


III. HYBRIDIZATION WITH A MICROWAVE 
CAVITY MODE 


C. Magnon linewidth 

Here we only consider magnons in insulating ferro- 
magnets. Metallic ones suffer from strong damping of 
magnons due to scattering by the conduction electrons. 

A number of magnon relaxation mechanisms are 
known in ferromagnetic insulators [24] [25] . At high tem- 
peratures around room temperature, magnon-magnon 
and magnon-phonon inelastic scatterings are domi- 
nant because of the large number of thermally excited 
magnons and phonons. The intrinsic magnon-magnon 
scattering is caused by the nonlinearity in the Holstein- 
Primakoff transformation [Eqs. (©- (0], while the latter 
is caused by the spin-lattice coupling [26] . However, both 
mechanisms are negligible at low temperatures we are in- 
terested in. 

At lower temperatures, extrinsic relaxation mecha¬ 
nisms become dominant. They are induced by the defects 
and impurities inside the crystals as well as the macro- 
scopic pores and roughness at the surfaces. In the inter- 
mediate temperature range, typically between 10-100 K, 
the linewidth often shows a large peak in the tempera¬ 
ture dependence. The peak height strongly depends on 
the amount of the defects such as rare-earth impurities 
and oxygen vacancies. It is known that the so-called slow- 
relaxation mechanism caused by the magnetic impurities 


In this section, we consider a sphere of ferromagnetic 
insulator embedded in a microwave cavity resonator. Re- 
lated experiments have been reported recently by a few 
other groups [3QH3^ . We study the coupling between the 
Kittel mode and a single discretized microwave mode in 
a cavity in the quantum limit. 

A. Theory 

The coupling between linear-polarized microwave pho- 
tons and the spin ensemble via the Zeeman effect is de- 
scribed by the Hamiltonian: 

Bint = 5MB F ■ ^o(ri) (a + a^) 

i 

= 5r/iBV^y]y]s„(ri) ■ Bo(rj) (a + a^), 

i n 

(17) 

where Bo(r^) is the linear-polarized microwave magnetic 
field in the cavity mode at the single photon level at the 
position of each spin. In the second line, we replace 
the Heisenberg operator with the sum of magnon op- 
erators multiplied by their orthonormal mode functions 
Sn(r^). Here, n is an index of the modes. 
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FIG. 3. Microwave cavity and YIG sphere. a) Photograph of a half of the microwave rectangular cavity made of oxygen-free 
copper. The cavity has dimensions of 22x18x3 mm and has the fundamental TEioi-mode frequency of 10.565 GHz. An YIG 
sphere is mounted in the cavity with a support rod made of aluminum oxide. The sphere is glued to the rod oriented to the 
crystal axis (110). b) Numerical simulation of the microwave magnetic-held distribution of the TEioi mode. The YIG sphere 
is located at the held maximum. 


We then replace the sum over the spins with a volume 
integral and apply the rotating wave approximation, fi- 
nally obtaining 


iiint = / drsn{r) ■ Bo(r) (a^c„ + dcl) , 

n 

(18) 

where V is the sample volume. For a cavity field spatially 
uniform within the sphere, we see from symmetry that 
the only mode with a finite coupling strength is the Kittel 
mode which has a spatially uniform function. In this case 
we obtain for Bq T z 


T^int — “1“ 


^eff = 




2h 


(19) 

( 20 ) 


where c = {l/y/‘^sN) is the annihilat ion operators 

of the Kittel mode, and Bq = |Bo(r)|. 


B. Yttrium iron garnet (YIG) 

In the following experiments, we use a single crystalline 
sphere of yttrium iron garnet (Y 3 Fe 50 i 2 ; YIG) as a fer- 
romagnetic sample. YIG is a celebrated ferromagnetic 
insulator [33] , used for various microwave devices includ- 
ing filters and oscillators. The absence of conduction elec- 
trons leads to its small spin-wave relaxation rate, which 
also makes YIG very attractive in spintronics applica- 
tions [34H36] . Strictly speaking, YIG is a ferrimagnetic 
material, but all the spins in a unit cell practically pre- 
cess in phase in the low energy limit, enabling us to treat 
it as ferromagnet. The net spin density 2sN/V in YIG 
is 2.1 X 10^^/iBcm“^, orders of magnitude higher than 
the typical numbers, 10 ^^ — IO^^/^b cm“^, in the param- 
agnetic spin ensembles used in quantum memory exper¬ 
iments. Thus, we expect strong interaction of the spin 
excitations with an electromagnetic field. 


C. Experiment 


To accommodate an YIG sphere in the confined space, 
where only a single electromagnetic mode exists in a 
certain frequency range, we use a three-dimensional mi¬ 
crowave cavity shown in Fig.[^. The picture shows a half 
cut of the cavity, and two pieces of them make a hollow 
cavity. The microwave magnetic-field distribution of the 
fundamental mode (rectangular TEioi) is simulated by 
COMSOL Multiphysics® (Fig. HD- A 0.5-mm-diameter 
YIG sphere is placed near the maximum of the magnetic 
field in order for obtaining the largest coupling strength 
and the uniformity of the field. 

We apply a static magnetic field of around 0.3 T by 
using a pair of neodymium permanent magnets and a 
10^-turn superconducting coii. They are connected in se¬ 
ries using a magnetic yoke made of pure iron (Japanese 
Industrial Standard SUY-1). The static field is oriented 
to the (100) crystal axis of the YIG sphere. The cavity 
has two SMA connectors for transmission spectroscopy. 
The center pins of the connectors are protruding into the 
cavity, such that their coupling strengths, and 

^out/27r, are about 0.5 MHz. We use a weak probe mi¬ 
crowave power of —123 dBm, which corresponds to the 
photon occupancy of 0.9 in the cavity mode. All the mea- 
surements are done in a dilution refrigerator: the ambi¬ 
ent temperature at the sample st age is 10 mK and the 
thermal photon/magnon occupancy at around 10 GHz is 
negligible. 

We measure the transmission coefiicient 5'2i(<^) of the 
cavity as a function of the probe frequency and the static 
magnetic field tuned by the bias current / in the coii 
(Fig. |§). A large normal-mode splitting is observed, 
manifesting strong coupling between the Kittel mode and 
the cavity mode. Cross sections of the intensity plot are 
shown in Fig. Et- At the degeneracy point where the Kit¬ 
tel mode frequency coincides with the cavity frequency, 
we see two largely-separated peaks. The peaks indicate 
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FIG. 4. Normal-mode splitting between the Kittel mode and the cavity mode TEioi. a) Amplitude of the transmission Re(S' 2 i) 
through the cavity as a function of the probe microwave frequency and the static magnetic field represented in the current I 
through the superconducting coii. The current I is defined to be zero at the anticrossing. The horizontal and diagonal dashed 
lines show the TEioi-mode and the Kittel-mode frequencies both obtained from the fitting. b) Cross sections at the static 
magnetic field corresponding to I — —3.5, —2.3, —1.1, and 0 mA. Solid curves are experimental data, and the dashed white 
lines are fitting curves based on the input-output theory. Eor clarity, the individual curves are offset vertically by 0.4 from 
bottom to top. c) Linewidth of the Kittel mode as a function of the temperature. The red dots show the linewidth obtained 
by fitting *S' 2 i(ct;) measured at each temperature. The dashed line is the fitting curve to the temperature dependence below 
1 K. The dashed bine line depicts the temperature-independent component of the fit. (Reprinted figure with permission from 
Y. Tabuchi, et al, Phys. Rev. Lett. 113, 083603 (2014). Copyright (2014) by the American Physical Society.) 


the hybridization of a Kittel-mode magnon and a cavity 
photon, i.e., formation of “magnon-polaritons”. Their 
decay rates are apparently much smaller than the cou- 
pling energy. 

We quantify the coupling strength the cavity-mode 
decay rates + /^out + ^int i and the Kittel-mode decay 
rate 7 ^, based on the model Hamiltonian in Eq. (20). 
Here, /^in, /^outi and /^int are the cavity decay rates 
through the input and output ports and the internal loss 
channel, respectively. Using the input-output theory, we 
derive the transmission coefiicient S 2 i{oj) as: 


In designing coupling strengths for various applica- 
tions, it is worth estimating the strength with numerical 
simulation. Figure shows the magnetic field distribu- 
tion Bq at the single photon level of the TEioi mode. The 
coupling strength can be easily calculated by the relation 
go = 7 H 0 / 2 . The simulated value of Hq = 5.5 pT/photon 
at the sample gives ^o/27r = 38.5 mHz, which agrees well 
with the experiment. 


E. Magnon linewidth 


521 M = 


y/ ^in^out 


,7/1 _ /1 i _ /^in+/^out+/^int I \9m \ 

LUc) 2 


( 21 ) 


The fitting curves, shown as the white dashed lines in 
Fig. [^, well reproduce the experimental data. From 
the fitting, we obtain = 47 MHz, (/^in + /^out + 

/^int)/27r = 2.7 MHz, and 7m/27r = 1.1 MHz. We find our 
magnon-cavity hybrid system deep in the strong coupling 
regime, ^ 7 , even in the quantum limit where the 

average photon/magnon number is less than one. 


D. Coupling strength 


The obtained coupling strength 5'm/27r of 47 MHz re- 
sults from the \/]V-enhancement according to Eq. (20). 


Given that the 0.5-mm-diameter sphere contains 1.4 x 
10 ^^ net spins, the single spin coupling strength is esti- 
mated to be 40 mHz. 


Little has been known about the linewidth of the Kittel 
mode in the temperature range attainable in a dilution 
refrigerator. We measure the temperature dependence of 
the resonance linewidth below 1 K (Fig. [^) and observe 
a peculiar behavior below 1 K; the linewidth is broadened 
as temperature decreases. 

The fitting curve based on the TLS model, as indi- 
cated with the green dashed line in Fig. [^, well agrees 
with the experimental data below 1 K. Note that the 
Kittel-mode frequency is used as a fixed parame- 
ter in the temperature-dependent term proportional to 
tanh(fkJni/ 2 /cBT). We also assume the presence of a 
temperature-independent contribution in the fitting. The 
parameters obtained imply that among the linewidth at 
the lowest temperature a fraction of 0.63 MHz is at- 
tributed to TLSs, and the other 0.39 MHz to surface 
scattering. The microscopic origin of the TLSs remains 
to be understood. The linewidth broadening above 1 K 
is ascribed to the slow-relaxation mechanism caused by 
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the magnetic impurities EU- 

An additional signature of the effect of the TLSs is 
found in the power dependence of the linewidth. Strong 
microwave drive saturates the TLSs coupled to the Kittel 
mode, resulting in the narrowing of the linewidth. Sim- 
ilar phenomena have been observed in superconducting 
microwave resonators interacting with TLSs [371 [38]. We 
confirmed experimentally that the Kittel-mode linewidth 
indeed became narrower at higher power of the drive. 
The power level causing the narrowing should be related 
with the dipole strengths and the relaxation rates of the 
TLSs. Comprehensive and systematic analyses of the 
linewidth with respect to the diameter and the quality 
of the spheres, the crystal orientations to the static field, 
and the Kittel mode frequency are also awaited. 


IV. COUPLING WITH A SUPERCONDUCTING 
QUBIT 

In the previous section we demonstrated coherent cou- 
pling of magnons to the cavity mode and investigated 
the linewidth in the quantum limit. However, that was a 
linear coupling between two harmonic oscillator modes. 
Although we observed the normal-mode splitting at the 
quantum limit near the ground state, there was no sig¬ 
nificant difference from the one we see in the classical 
limit, e.g., at room temperature and with much stronger 
probe microwave. The correspondence principle makes it 
difiicult to distinguish quantum and classical behaviors 
in such a linear system. Nonlinearity, or aharmonicity, 
is required for a ciear demonstration of the quantum be- 
havior. Thus, in this section, we move one step further 
by introducing a well-controllable two-level system, i.e., 
a superconducting qubit. 


A. Theory 

In order to implement a coupling between a supercon¬ 
ducting qubit and a magnon in YIG, we exploit the cavity 
quantum electrodynamics architecture. In the scheme il- 
lustrated below, two heterogeneous systems — the qubit 


and the Kittel mode — are linked through electromag- 
netic fields in a microwave cavity resonator. 

We use an elongated cavity shown in Fig.[^, which si- 
multaneously accommodates a qubit and an YIG sphere. 
The 5-cm-long cavity has eigenmodes of TEiop {p = 
1,2, •••) as the lowest-frequency modes. The resonant 
frequencies ujiop are determined by the width W and 
length L of the cavity, and denoted as: 



where cq is the speed of light. In Fig. It the simulated 
electric and magnetic field distributions of the TE 102 
mode are shown. 

We use a transmon qubit which has two aluminum 
pads and a Josephson junction bridging them. The qubit 
chip is embedded in the cavity as seen in the bottom pic¬ 
ture of Fig. It- The submillimeter-sized dipole antenna 
consisting of the two pads electrically couples to the mi¬ 
crowave field in the cavity. The qubit is placed where 
the electric field is close to the maximum. On the other 
hand, the Kittel mode of the YIG sphere magnetically 
couples to the same cavity mode. The corresponding 
system Hamiltonian is written as: 


y-sys/fi = Y 


+ OJm C^C, 


(23) 


where uJiop is the resonant frequency of TEiop, cjq and 
cjrn are the qubit and the Kittel mode frequencies, and 
dp, b = + 110qq(^ + M ^ annihilation 

operators for the cavity mode, the qubit, and the Kittel 
mode, respectively. The transmon qubit is represented 
as an anharmonic oscillator with an anharmonicity of 
(a < 0. The lowering operator for the qubit subspace 
is defined as a~ = [cr+]’*' = Vh^ where V = |0)(0| -1- 
|1)(1| is a projection operator. The annihilation operator 
of the Kittel mode c = (l/\/2sY) stems from 

the Holstein-Primakoff transformation. The interaction 
Hamiltonian is written as: 





^ ^ 9 m,i 0 p -|- h.c. 


(24) 


where g^pop and gmpop are the coupling strengths of the 
cavity TEiop mode to the qubit and the Kittel mode, re¬ 
spectively. Note that the direct interaction between the 
qubit and the Kittel mode is negligible. In the large 
detuning regime where ^m,i 0 p < |cjq - < 


\uji 0 p — ^q|, l^iop — ^m|, the interaction Hamiltonian can 
be rewritten in the rotating frame of the system Hamil¬ 
tonian under a perturbative treatment up to the first 
order [39] , and expressed as: 
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a) 




FIG. 5. Qubit-cavity-YIG hybrid System, a) Photographs of an YIG sphere and a superconducting qubit mounted in an 
extended cavity resonator. The YIG sphere is mounted in one side of the cavity. On the other side, the superconducting 
transmon qubit fabricated on a Silicon chip is seen in the deep cavity trench. The cavity with dimensions of 25 x 3 x 53 mm 
supports the lowest-frequency (TEioi, TE 102 , and TE 103 ) modes with the frequencies of 6.987 GHz, 8.488 GHz, and 10.461 GHz, 
respectively. b) Numerical simulation of the TE 102 mode. Intensities of the microwave electric and magnetic fields are shown, 
respectively; the bine arrows indicate the directions of the helds. The qubit chip is placed near the electric held antinode. 
The qubit with a millimeter-sized dipole antenna strongly couples with the electric held. The YIG crystal is placed where the 
magnetic held is large. The color scales indicate held intensities at the single photon level. c) Schematic energy diagram of the 
qubit-cavity-YIG hybrid System. The left shows unperturbed (bare) frequencies of the qubit, the cavity and the Kittel modes. 
The cavity and the Kittel modes are subject to frequency shifts due to the coupling with the qubit, as indicated in the middle 
of the diagram. The cavity modes induce the Lamb shift of the qubit, as depicted on the right-hand side. 


nj /fc 9m.,10p 

Ttint/^ = 2^ - CiLci' 


9in,10p /sf r 
-1- 1. - I r 


. LLp ^ ^ I 0 qq(^l 


+ Y (Xl0p|0)qq(0| + YXlOpl0qq(^lj «lop^lOp “ Y ^ ^Op^lOp C, (25) 


Z>0 


where xiop = \ 9 q,i 0 p\‘^/{^lOp - ^q), >^iop = 

^ |^q,10p| /[^lOp ~ ^q ~ ~ XlOp ~ |^q,10p| [(^ Y 

l)/(<^i 0 p~^q—^ — — <^)]- The hrst and 

the third terms are dispersive shifts due to the coupling 
between the cavity modes and the Kittel mode. The sec- 
ond term indicates the Lamb shift of the qubit frequency 
arising from the coupling to the cavity mode. The Lamb 
shift for the first excited States coincides with XiOp- 
The fourth term shows the qubit-state-dependent cavity 
frequency shift or the photon-number-dependent qubit 
frequency shift. It is worth noting that the last term in¬ 
dicates the static interaction between the cavity modes 
and the Kittel mode, which originates from the fact that 
the spin ensemble is not a perfect bosonic system. Be- 
cause of the factor l/N, however, for large N such effect is 
not observable with usual experimental parameters. The 
major energy shifts are sum marized in Fig. [5]:. In the 
regime where gq,iop9m,iOp /Vkiop - Wql |wiop - ui^l < 


|wq - Wnil < |wiop - Wql, \LOiop - Wni | , there is a qubit- 
state-dependent frequency shift ^ of the Kittel mode m- 
Although it appears only in the third-order perturbative 
treatment, the shift is stili observable when the Kittel- 
mode and the qubit frequencies are close enough to meet 
the frequency condition. Such coupling is usable for 
counting the magnon number in the Kittel mode via a 
Ramsey interferometry using the qubit, for example. 


The coupling between the qubit and the Kittel mode 
is mediated by the cavity modes TEiop when the qubit 
and the Kittel-mode frequencies are degenerate with each 
other and detuned from the cavity so that jcjq — 

9 q, 10 p'> 9 m, 10 p YC l^lOp ^q| — l^lOp | ‘ 

rewrite the interaction Hamiltonian in Eq. (24) in the 
corresponding rotating frame by adiabatically eliminat- 
ing the cavity modes: 


'Hint/fi = 9q-m O + h.C., 


(26) 





























where 5q_ni = Ep 5m,i0p 5q,i0p/(‘^i0p -‘^q)- It is in- 
terpreted as that the qubit and the Kittel mode ex- 
change their energy quanta, i.e., the qubit excitation and 
a magnon, through the virtual excitation of the cavity 
mode TEiop. 


B. Experiment 


branches. When we apply a resonant microwave to the 
branches, they are excited and detected as a phase shift 
in the reflection coefficient of the cavity probe microwave 
at the TEios frequency, given that each branch has a 
fraction of the qubit wave function. In this scheme which 
keeps the cavity mode TE 102 empty, unwanted dephasing 
caused by photon number fluctuations in the “coupler” 
mode can be avoided. 


Here we focus on the roles played by the cavity modes 
TE 102 and TE 103 (Eig. [^). The cavity mode TE 102 is 
utilized for mediating the virtual-photon interaction be- 
tween the qubit and the Kittel mode, while the cavity 
mode TE 103 is used for the qubit measurement. Al- 
though the Kittel mode also couples to the TE 103 mode, 
we can perform a selective qubit measurement via the 
mode because the magnon-number-dependent shift of the 
last term in Eq. (25) is sufficiently small and safely ig- 
nored. We introduce probe and excitation microwaves 
into the cavity through an SMA port. The external cou- 
pling strengths and internal losses of each cavity mode are 
^los/^TT = 2.75 MHz, /^io3,int/27r = 1.26 MHz, /i:io2/27r 
= 0.55 MHz, and /^io3,int/27r = 1.73 MHz, respectively. 

A 0.5-mm-diameter YIG sphere is placed where the 
microwave magnetic field of the cavity mode TE 102 is 
substantial, as indicated in Eig. It- The field is paral- 
lel to the (110) crystal axis of the YIG sphere. We also 
apply a static magnetic field locally to the sphere in par- 
allel with its (100) crystal axis. A doubly-layered mag¬ 
netic shield protects a superconducting qubit from the 
stray magnetic field which is approximately a few tens 
of gauss at the qubit position. The coupling strength 
^m,i 02 of the Kittel mode to the cavity mode TE 102 
and the decay rate 7 ^ of the Kittel mode are obtained 
from the anticrossing spectrum in the reflection mea¬ 
surement: po2/27r = 21 MHz and 7ni/27r=1.8 MHz. 

The transmon qubit has the dressed resonant frequency 
a;q/27r = (cjq — X]qXi0p)/27r of 8.158 GHz and the an- 
harmonicity (a/27r of —158 MHz. A separate time-domain 
measurement shows that the energy relaxation time Ti = 
273 ns of the qubit is almost at the Purcell limit where 
l/J^i = Ep XiOp(kiop + Ki0p,i„t)/(wi0p - Wq). We deduce 

the Lamb shift = Xiop from the frequency differ- 

ence between the bare- and the dressed-cavity frequen- 
cies, which is accurate up to the first-order perturbative 
treatment; they are evaluated to be Xio2/27r = 75 MHz 
and Xio3/27r = 9 MHz. The coupling strengths of the 
qubit to the cavity modes are 5'q,io2/27r = 117 MHz and 
^q,i03/27r = 141 MHz. The qubit-state-dependent cavity 

shift (Xio 3 “ Xi03)/27r of 1.2 MHz is used for the qubit 
spectroscopy. The qubit linewidth 7q/27r is 2.0 MHz. 

Eigure illustrates the measurement scheme of the 
hybridized qubit-magnon modes. At the degeneracy 
point where the Kittel-mode frequency coincides with the 
qubit frequency, the virtual excitation of the cavity mode 
TE 102 mediates interaction between them as expressed 
in Eq. (26), so that the qubit frequency splits into two 


C. Resuit and discussion 


Eigure shows the reflection coefhcient at the TE 103 
frequency as a function of the excitation frequency and 
the coii current. The pronounced splitting of the qubit 
frequency at the degeneracy point manifests magnon- 
vacuum-induced Rabi splitting, which indicates coherent 
coupling between the qubit and the Kittel mode. The ob¬ 
tained coupling strength ^q_m/27r is 11.4 MHz, which is 
much larger than the qubit decay rate 7 q and the Kittel- 
mode decay rate 7 m- Thus the hybridized qubit-magnon 
System also stands deep in the strong coupling regime 
where ^q_m > 7q^7m- 

By using Eq. (26), we estimate the coupling strength 


for a detuning W 102 - Wq = W 102 - Wq - Ep=i,2,3 
= 183 MHz and obtain 13.4 MHz. The discrepancy of 
17% between the calculated and the measured strengths 
is ascribed to the first-order approximation of the Hamil- 
tonian; the value |^q,io 2 /(<^i 02 “^q)| ^ 0 -b should be far 
less than unity for a convergence of the series. Another 
possible factor is the estimation error in the bare qubit 
frequency cjq. Because all the cavity modes contribute 
to the Lamb shift of the qubit, it is not straightforward 
to accurately estimate the shift in a multi-mode cavity, 
especially when the modes are crowded. 


V. CONCLUSIONS AND OUTLOOK 

We demonstrated coherent coupling between a single 
magnon excitation in a ferromagnet and a microwave 
photon in a cavity as well as a superconducting qubit. 
It was proven that the uniform spin precession mode, 
i.e., the Kittel mode, of a millimeter-scale ferromagnetic 
sphere behaves quantum mechanically: the qubit showed 
a Rabi splitting induced by the vacuum fluctuations of 
the Kittel mode. 

The technique developed here exploits the advanced 
Circuit-QED and microwave quantum optics technologies 
based on superconducting qubits. It enables generation 
and characterization of non-classical States of magnons 
and thus opens a new field of quantum magnonics. It 
is also considered to be the ultimate limit of spintronics. 
We expect that further studies will reveal the dissipation 
mechanisms in the single magnon regime. 

The coherent transfer of the qubit state to the magnon 
mode suggests a possible link to quantum Information 
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FIG. 6. Magnon-qubit hybridization. a) Energy level diagram illustrating the coupling scheme. The superconducting transmon 
qubit is detuned from the TE 102 cavity mode by —178 MHz, and the Kittel mode is almost degenerate with the qubit. The 
couplings of the TE 102 mode with the qubit and the Kittel mode induce their interaction via the virtual-photon excitation; the 
qubit-magnon degeneracy is lifted by 2^q_m- The TE 103 cavity mode is dispersively coupled to the qubit so that its resonant 
frequency depends on the qubit state, b) Magnon-vacuum-induced Rabi splitting. The change in the real part of the reflection 
coefficient at the TE 103 frequency, Re(Ar), is plotted as a function of the qubit excitation frequency and the static magnetic 
field which is locally applied to the YIG sphere. The static field is represented in the coii current /, which is defined to be zero 
at the anticrossing. The white solid curve shows the cross section at / = 0 mA. The powers used for probing the TE 103 mode 
and for exciting the qubit and the Kittel mode are —135 dBm and —130 dBm, respectively. Note that the measurement here 
was done with the same sample set as in Ref. [19]. By readjusting the alignment of the YIG sphere, we were able to remove 
the effect of other magnetostatic modes previously observed. 


networks in optics. In contrast to superconducting Cir¬ 
cuit s, spins in insulating crystals can interact coherently 
with light, as demonstrated in quantum memory exper- 
iments in the optical domain [411 IU]- In those experi- 
ments, paramagnetic spin ensembles have been used. It 
is of great interest to look into coherent coupling between 
magnons in ferromagnet and light. 
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